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I. INTRODUCTION 



Hadronic matter under extreme conditions has attracted a lot of interest in recent years. 
On the one hand, many theoretical and experimental efforts have been devoted to the 
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discussion of heavy ion collisions where the temperature is high. On the other hand, the 
physics of neutron stars has become a hot topic which connects astrophysics with high 
density nuclear physics. In 1934, Baade and Zwicky [1] suggested that neutron stars could 
be formed in supernovae. The first theoretical calculation of a neutron star was performed 
by Oppenheimer and Volkoff [2], and independently by Tolman [3]. Observing a range of 
masses and radii of neutron stars will reveal the equations of state (EOS) of dense hadronic 
matter. Determination of the EOS of neutron stars has been an important goal for more 
than two decades. Six double neutron-star binaries are known so far, and all of them have 
masses in the surprisingly narrow range 1.36±0.08M sun [4,5]. A number of early theoretical 
investigations on neutron stars were based on the non-relativistic Skyrme framework [6]. 
Since the Walecka model [7] was proposed and applied to study the properties of nuclear 
matter, the relativistic mean field approach has been widely used in the determination of 
the masses and radii of neutron stars. These models lead to different predictions for neutron 
star masses and radii [8,9]. For a recent review see Ref. [10]. Though models with maximum 
neutron star masses considerably smaller than 1.4M sun are simply ruled out, the constraint 
on EOS of nuclear matter (for example, the density dependence of pressure of hadronic 
system) has certainly not been established from the existing observations. 

In the process of neutron star formation, /^-equilibrium can be achieved. As a conse- 
quence, hyperons will exist in neutron stars, especially in stars with high baryon density. 
These hyperons will affect the EOS of hadronic matter. As a result, the mass-radius rela- 
tionship of strange hadronic stars will be quite different from that of pure neutron stars. 
The simplest way to discuss the effects of hyperons is to study strange hadronic stars in- 
cluding only A hyperons. This is due to the fact that A is the lightest hyperon and the A-N 
interaction is known better than other hyperon-nucleon interactions. However, one must 
also consider hyperons with negative charge in neutron stars because the negatively charged 
hyperons can substitute for electrons. There have been many discussions of strange hadronic 
stars including A hyperons, A and S~ or even the whole baryon octet [11] - [20]. 

At high baryon density, the overlap effects of baryons are very important and the quark 
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degrees of freedom within baryons should be considered. There are some phenomenological 
models based on the quark degrees of freedom, such as the quark meson coupling model [21], 
the cloudy bag model [22], the quark mean field model [23] and the NJL model [24]. Several 
years ago, a chiral SU(3) quark mean field model was proposed [25,26]. In this model, 
quarks are confined within baryons by an effective potential. The quark-meson interaction 
and meson self-interaction are based on SU(3) chiral symmetry. Through the mechanism of 
spontaneous symmetry breaking the resulting constituent quarks and mesons (except for the 
pseudoscalars) obtain masses. The introduction of an explicit symmetry breaking term in 
the meson self-interaction generates the masses of the pseudoscalar mesons which satisfy the 
partially conserved axial-vector current (PCAC) relations. The explicit symmetry breaking 
term in the quark-meson interaction gives reasonable hyperon potentials in hadronic matter. 
This chiral SU(3) quark mean field model has been applied to investigate nuclear matter 
[27], strange hadronic matter [25], finite nuclei, hypernuclei [26], and quark matter [28]. 
Recently, we improved the chiral SU (3) quark mean field model by using the linear definition 
of effective baryon mass [29]. This new treatment is applied to study the liquid-gas phase 
transition of asymmetric nuclear system and strange hadronic matter [30,31]. By and large 
the results are in reasonable agreement with existing experimental data. 

In this paper, we will study the neutron star and strange star in the chiral SU (3) quark 
mean field model. The paper is organized in the following way. In section II, we briefly 
introduce the model. In section III, we apply this model to investigate the neutron star and 
strange hadronic star. The numerical results are discussed in section IV. We summarize the 
main results in section V. 

II. THE MODEL 

Our considerations are based on the chiral 577(3) quark mean field model (for details 
see Refs. [25,26]), which contains quarks and mesons as the basic degrees of freedom. In 
the chiral limit, the quark field ^ can be split into left and right-handed parts ty L and fy R : 
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# = m L + m R . Under SU(3) L x SU(3) R they transform as 



R 



(1) 



The spin-0 mesons are written in the compact form 
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M 
Mt 



(2) 



s±Ti = ^=]T( s a ±^)A a , 

where s a and p a are the nonets of scalar and pseudoscalar mesons, respectively, \ a (a = 
1, 8) are the Gell-Mann matrices, and A = ^/| /. The alternatives indicated by the plus 
and minus signs correspond to M and Aft, respectively. Under chiral SU(3) transformations, 
M and Aft transform a. M — W LME) and Aft -> Aft' = i^M+Lt. The spin-1 mesons 
are arranged in a similar way as 



1 1 8 

M = - (Va ± A„) = — = y (v a u ± A c 



(3) 



with the transformation properties: Z M — > = Ll^lJ , r M — > = Rr^B) . The matrices E, 
II, U M and can be written in a form where the physical states are explicit. For the scalar 
and vector nonets, we have the expressions 

\ 
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s = —j= y s a x a 

v2 a=0 



if* 



*0 



K* 



K 



*o 



c 



(4) 



p-, ^{^ -pi) K 



(5) 



Pseudoscalar and pseudovector nonet mesons can be written in a similar fashion. 
The total effective Lagrangian is written: 



(6) 



where £o = i^fj^d^ is the free part for massless quarks. The quark-meson interaction C q M 
can be written in a chiral SU (3) invariant way as 



= ^jf (E *A + *7 5 X>«A.) * - ^ (V E <A a - 7 V E«^) *■ (7) 

From the quark-meson interaction, the coupling constants between scalar mesons, vector 
mesons and quarks have the following relations: 

9 s _ n u _ d _ u _ d _ _ |_ s s _ s _ u _ d _ n /o\ 

- 9a - ~9a - 9a - 9a - ■ ■ ■ - 7y|#<> 9a ~ 9a ~ 9c, ~ 9( ~ U , (8 J 

9v __ u d u „d ^ „s „s „s „u „d 



, )sJ1 -9 u P = -g d P = 9l = gt = ... = £ = 91 = g$ = 4 = o. (9) 

In the mean field approximation, the chiral-invariant scalar meson £ ss and vector meson 
Lyy self-interaction terms are written as [25,26] 

£ee = -\ hx 2 (- 2 + C 2 ) + h (V + C 2 ) 2 + h + C 4 ) + h X a\ 



-h X 4 - TX 4 ln^ + ^X 4 ln^>, (10) 
4 Xo 3 a Co 

£vv = \ 4 K^ 2 + ™> 2 + ™5^) + 94 U + 6c 2 p 2 + P 4 + 20 4 ) , (11) 
^ Xo 

where 5 = 6/33; <t and ( are the vacuum expectation values of the corresponding mean 
fields a, ( which are expressed as 

<r = -F«, (o = ^(F w -2F K ). (12) 

The vacuum value xo is about 280 MeV in our numerical calculation. The Lagrangian £ x sb 
generates the nonvanishing masses of pseudoscalar mesons 

Xo 

leading to a nonvanishing divergence of the axial currents which in turn satisfy the partial 
conserved axial-vector current (PCAC) relations for it and K mesons. Pseudoscalar, scalar 
mesons and also the dilaton field x obtain mass terms by spontaneous breaking of chiral 
symmetry in the Lagrangian of Eq. (10). The masses of u, d and s quarks are generated 
by the vacuum expectation values of the two scalar mesons a and (. To obtain the correct 
constituent mass of the strange quark, an additional mass term has to be added: 



s* X 

L XSB = ~2 



mlF^a + (V2m 2 K F K - ^F^j ( 



(13) 



£Am a = -Am s qSq (14) 

where S — | (i — XsV^j = diag(0, 0, 1) is the strangeness quark matrix. Based on these 
mechanisms, the quark constituent masses are finally given by 

m u = m d = --^=<7 and m s = -g s ( + Am s . (15) 

The parameters g s = 4.76 and Am s = 29 MeV are chosen to yield the constituent quark 
masses m q = 313 MeV and m s = 490 MeV. In order to obtain reasonable hyperon potentials 
in hadronic matter, we include an additional coupling between strange quarks and the scalar 
mesons a and ( [25]. This term is expressed as 

C h = [h 1 (a-<j ) + h 2 ((-( )]ss. (16) 

Therefore, the strange quark scalar-coupling constants are modified and do not exactly 
satisfy Eq. (8). The hyperon potentials were listed in our previous paper [29]. In the quark 
mean field model, quarks are confined in baryons by the Lagrangian C c = — ^ Xc ^ (with 
Xc given in Eq. (17), below). We note that this confining term is not chiral invariant. 
Possible extensions of the model which would restore chiral symmetry in this term have 
been discussed in Ref. [32]. 

The Dirac equation for a quark field ^ under the additional influence of the meson 
mean fields is given by 

-ia ■ V + PXc(r) + pm* 

where a = 7°7 , (3 = 7 , the subscripts % and j denote the quark % {% = u, d, s) in a baryon of 
type j (j = N, A, E, S) ; xdj') is a confinement potential, i.e. a static potential providing the 
confinement of quarks by meson mean-field configurations. In the numerical calculations, we 
choose Xc(r) — \k c r 2 , where k c — 1 (GeV fm~ 2 ), which yields baryon radii (in the absence 
of the pion cloud [33]) around 0.6 fm. The quark mass m* and energy e* are defined as 

m* = -gi<r - g\C, + m i0 (18) 
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(17) 



and 



e* = et- g^-g l p p -g\4>, (19) 

where is the energy of the quark under the influence of the meson mean fields. Here 
m i0 — for i — u, d (nonstrange quark) and m i0 = Am s for i = s (strange quark). The 
effective baryon mass can be written as 

M* = J2n ije *-El (20) 

i 

where riij is the number of quarks with flavor "i" in a baryon with flavor j, with j = 
N {p, n} , E {E^, S } , S S - } , A and E® was found to be only very weakly dependent 
on the external field strength. We therefore use Eq. (20), with E® a constant, independent 
of the density, which is adjusted to give a best fit to the free baryon masses. Compared with 
the earlier square root ansatz, here we use the linear definition of effective baryon mass. 
As we have explained in Ref. [29] the linear definition of effective baryon mass has been 
derived using a symmetric relativistic approach [34], while to the best of our knowledge, no 
equivalent derivation exists for the square root case. 



III. HADRONIC SYSTEM 

Based on the previously defined quark mean field model the thermodynamical potential 
for the study of hadronic systems is written as 

n= V /" d3fc fin (l + e" W>"^^ 

J= W,H ( 27r ) 3 J ° 1 V 1 V l] 

(21) 

where E*(k) = M* 2 + k 2 and M* is the effective baryon mass. The quantity Uj is related 
to the usual chemical potential /ij by Vj = fij — g^u — g 3 p p — gfyj). The mesonic Lagrangian 

£m = ^ss + £-vv + £-xSB (22) 
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describes the interaction between mesons which includes the scalar meson self-interaction 
£ee, the vector meson self-interaction C V v an d the explicit chiral symmetry breaking term 
£-xSB defined previously in Eqs. (10), (11) and (13). The Lagrangian Cm involves scalar 
(a, ( and x) an d vector (cu, p and fa) mesons. The interactions between quarks and scalar 
mesons result in the effective baryon masses M*. The interactions between quarks and 
vector mesons generate the baryon-vector meson interaction terms. The energy per volume 
and the pressure of the system can be derived as e = f2 — ^§p + VjPj and p = —Q, where pj 
is the density of baryon j. At zero temperature, fl can be expressed as 

1 



n = - £ 



j=N,A,x; 



24tt 2 





v) - Mf 


1/2 


2z/| - 5M* 2 " 











+ 3M* 4 ln 



1/2 



M* 



Cm, 



(23) 



The equations for mesons fa can be obtained by the formula j^p = 0. Therefore, the 
equations for a, ( and x are 

25 4 , X 2 



k oX 2 a - 4fci (a 2 + C 2 ) or - 2k 2 a 3 - 2k^C - 
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,Xo/ da 

k oX 2 ( - 4*! (a 2 + C 2 ) C - 4A; 2 C 3 - A^a 2 - ^X 4 + ^ (>/2ro^ fc - 



(24) 



,Xo, 



< ^j^j >= 0, 



(25) 



^oX (o" 2 + C 2 ) - h& 2 ( + ( 4A; 4 + 1 + 4/n^ - -^/n^ 
v 7 V Xo 3 (T 2 C 0/ 



+ 



2x 
Xo 



X / 2 2 

Xo 



2 2 

m p 



m 2 2 ) 



0, 



(26) 



where < ipjipj > is expressed as 



with k F . = yVj - Mf. 

For the /3-equilibrium, the chemical potentials for the baryons satisfy the following equa- 
tions: 

Pa = Ps° = Ps° = P« = P P + Pe = P P + p M , (28) 
// s + = (29) 

Ats- = Ps- = P« + /v (30) 

There are only two independent chemical potentials which are determined by the total 
baryon density and neutral charge: 

Pb = P P + pn + Pa + Ps+ + Ps« + Ps- + Ph» + Pa- , (31) 

P P + Ps+ - Ps- - Pb- - Pe - P M = 0. (32) 

In order to get the mass-radius relation, one has to resolve the Tolman-Oppenheimer- 
Volkoff (TOV) equation: 

dp _ \p(r) + e(r)] [M(r) + 4irr 3 p(r)] 



dr r(r - 2M(r)) 

where 



(33) 



Mir) =4tt f r e(r)r 2 dr. (34) 
Jo 

With the equations of state, the functions, such as M(r), p{r) andp(r), etc. can be obtained. 

IV. NUMERICAL RESULTS 

The parameters of this model were determined by the meson masses in vacuum and the 
saturation properties of nuclear matter which were listed in the table 1 of Ref. [29]. The 
improved linear definition of effective baryon mass is chosen in our numerical calculations. 



We first discuss the equations of state of neutron matter and strange hadronic matter which 
are needed for the calculation of neutron stars. For pure neutron stars, there are only 
neutrons present. For strange hadronic stars, with increasing baryon density, other kinds of 
baryons will appear. In Fig. 1 we show the fractions of octet baryons versus density with 
/9-equilibrium. With the increasing of baryon density, the neutron fraction decreases slowly 
from 1. If the density is lower than about 0.19 fm~ 3 , the fraction of electrons is the same 
as that of protons which makes the system charge neutral. The muon appears when the 
density is in the range 0.19 - 0.98 fuT 3 . The maximum fractions of muons and electrons 
appear at ps — 0.4 fm~ 3 . Their fractions decrease with the increasing fractions of hyperons. 
When the density is larger than about 0.4 fm~ 3 , the S~ hyperons appear and the fraction 
of neutrons decreases faster. After the density is larger than about 0.57 fm -3 , A hyperons 
start to appear. The fraction of S~ hyperons decreases with the increasing density after 2~ 
hyperons appear where the density is about 0.84 fm~ 3 . 
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FIG. 1. The fractions of proton, neutron, A, E and E of strange hadronic stars versus baryon 
density with /3-equilibrium. 
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FIG. 2. The effective baryon masses and meson mean fields versus baryon density with 
/^-equilibrium. 

The density dependence of the effective baryon masses and scalar mean fields are shown 
in Fig. 2. The a field decreases quickly with the increasing baryon density when the density 
is small, p B < 0.4 fm~ 3 . This is because at small baryon density, the nucleon is dominant 
and there are no hyperons. With the increasing of density, more and more hyperons appear. 
As a result, the ( field decreases quickly. At a broad range of densities, the value of x 
changes little. 

In Fig. 3, the pressure versus baryon density is shown. The dashed and solid lines are 
for the pure neutron star and the strange hadronic star with /9-equilibrium, respectively. 
When the density is low, the two curves are close to each other. With the increasing of 
baryon density, the contributions of protons and hyperons are not negligible. The inclusion 
of hyperons will soften the equation of state of hadronic matter. As a result, at a given 
baryon density the pressure of strange hadronic matter is smaller than the corresponding 
pressure of pure neutron matter. The pressure p versus energy density e is shown in Fig. 4. 
Again, one can see that the equation of state of strange hadronic matter is softer than that 
of pure neutron matter. 




P 8 ( fm ) 

FIG. 3. The pressure of hadronic matter p versus baryon density ps- The dashed and solid 
curves are for pure neutron stars and strange hadronic stars with /3-equilibrium, respectively. 
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FIG. 4. The pressure of hadronic matter p versus energy density e. The dashed and solid curves 
are for pure neutron stars and strange hadronic stars with /3-equilibrium, respectively. 



We now study neutron stars with the obtained EOS. By solving the TOV equation, the 
baryon density versus radius can be obtained which is shown in Fig. 5. The central densities 
p c are chosen to be 3po and 5po where po (0.16 fm~ 3 ) is the saturation density of symmetric 
nuclear matter. The dashed and solid lines are for pure neutron stars and strange hadronic 
stars with /3-equilibrium, respectively. With the increasing radius, the density of strange 
hadronic stars decreases a little faster than that of pure neutron stars which results in a 
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smaller radius. The radii of stars are not sensitive to their the central density. For example, 
for p c of 3p and 5p , the radii are both around 11-12 km. 



E 




R ( km ) 

FIG. 5. The baryon density of a hadronic star versus radius. The dashed and solid curves are 
for pure neutron stars and strange hadronic stars with /3-equilibrium, respectively. 

We plot the star mass ratio M/M sun versus central baryon density in Fig. 6. The max- 
imum mass of pure neutron stars is about 1.8M sun with a central density 1.05 fm~ 3 . After 
the central density is larger than 1.05 fm~ 3 , the star will become unstable. The maximum 
mass changes to 1.45M sun when hyperons are included. In the range 3p < Pc < ®Po, the 
masses of pure neutron stars and strange hadronic stars are 1.48M sun < M < 1.8M sun and 
1.23M sun < M < 1.45M sun , respectively. Our results are reasonable compared with the ob- 
servation of the six known stars with masses in the range 1.36±0.08M sun , since the "neutron 
star" is in fact a strange hadronic star with /^-equilibrium. We should also keep in mind that 
there are some heavy stars reported in recent years. For PSR J0437-4715, the mass is found 
to be 1.58 ± 0.18M sun [35]. For Vale X-l, Cygnus X-2 and 4U 1820-30, their masses are 
determined to be 1.87lg;??M sun [36], 1.8±0.4M sun [37] and ~ 2.3M sun [38,39]. The rotation 
of a star can increase its mass by ~ 10% [40]. Therefore, the calculated maximum mass of 
strange hadronic stars can be as large as 1.6M sun . If the heavy stars such as 4U 1820-30 are 
confirmed, the strange hadronic star would be ruled out if this model is a good description 
of Nature. It is possible to increase the maximum star mass by making the EOS stiffer at 

13 



higher densities. Whether the inclusion of a quark core in the strange star will result in a 
large-maximum mass is an interesting topic. 



2.0-r 

1.5- 

2 1.0- 
0.5- 
0.0- 

0.0 0.4 0.8 1.2 1.6 

P B ( f m 3 ) 

FIG. 6. The masses of hadronic stars versus their central baryon densities. The dashed and solid 
curves are for pure neutron stars and strange hadronic stars with /3-equilibrium, respectively. 

In Fig. 7, the masses of stars versus their radii are shown. For pure neutron stars, when 
their masses are in the range 0.5M sun < M < 1.8M sun , their radii are about 11.0-12.3 km. 
For the strange hadronic stars, when the masses are in the range 0.5M sun < M < 1.45M sun , 
the radii are about 10.7-11.7 km. With the same mass (M > 0.2M sun ), strange hadronic 
stars have smaller radii compared with pure neutron stars. Because the size of neutron stars 
is small, it is very difficult to observe and measure their radii directly. Different indirect 
methods lead to different values of radii with large errors. For example, for the RX J1856- 
3754, the radius varies from 5 km to 15 km with a mass of 1.4M sun [5]. More accurate values 
are needed to obtain a more strict constraint on the EOS of hadronic matter. 
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FIG. 7. The masses of hadronic stars versus their radii. The dashed and solid curves are for 
pure neutron stars and strange hadronic stars with /3-equilibrium, respectively. 



V. SUMMARY 

We have investigated pure neutron stars and strange hadronic stars in the chiral SU(3) 
quark mean field model. The A, E and S hyperons are included in the model. The proton 
and hyperon contributions to the system are important at high baryon density when (5- 
equilibrium is achieved, and soften the EOS of hadronic matter. The maximum pure neutron 
star mass is about M = 1.8M sun with a corresponding radius R = 11.0 km and central 
density p c = 1.05 fm~ 3 . For the strange hadronic stars, the maximum masses are about 
1.45M sun and the corresponding radii and central density are R = 10.9 km and p c = 1.0 
fm -3 . When the central densities are between 3po and 6po, the masses of stars are in the 
range 1.23M sun < M < 1.45M sun (strange hadronic stars) and 1.48M sun < M < 1.8M sun 
(pure neutron stars). If the masses of stars are larger than 0.5M sun , the typical values of 
radii are 10.7-11.7 km (strange hadronic stars) km and 11.0-12.3 km (pure neutron stars). 

Our results are reasonable compared with astrophysical observations where the six known 
neutron stars have masses in the narrow range 1.36 ± 0.08M sun . Accurate values of radii 
for neutron stars are needed to get a more strict constraint on the EOS of hadronic matter. 
As for the heavy stars, for example 4U 1820-30, if its mass M ~ 2.3M sun is confirmed, 

15 



then strange hadronic stars are obviously ruled out if the model explored herein is a good 
description of Nature. It is therefore of interest to see whether including quark degrees of 
freedom can lead to this large mass. 
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